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1 Introduction 

It has long been asserted that the number of ways to combine six 2 x 4 LEGO 
blocks of the same color is 

102981500 

This number was computed at LEGO in 1974 ([2]) and has been systematically 
repeated, for instance in 0] p. 15] and [Hj- Consequently, the number can be 
found in several "fun fact" books and on more than 250 pages on the World 
Wide Web 1 . However, this number only gives (with a small error, as we shall 
see) the number of ways to build a tower of LEGO blocks of height six. The 
total number of configurations is 

915103765 

as found by independent computer calculations by the second author and by 
Abrahamsen p^. This figure has now been accepted by LEGO Company, cf. 



oooooool oooon 



^^^^ 



looool loop 



Igggg g lg gggl 



joooo 



40 SSI si 



ooo oooo ooo|o o| 



Figure 1: 46 basic positions 

*LEGO is a trademark of LEGO Company 
1 Google search October 2004 
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We consider contiguous buildings of LEGO blocks, disregarding color, and 
identify them up to translation and rotation. We think of a 6 x w LEGO block 
as a subset of M 3 of the form 

[x; x + 6] x [y, y + w] x [z; z + 1] 

or 

[x; x + w) x [y, y + b] x [2; z + 1] 

Since the top and bottom of a LEGO block are distinguishable we will only 
consider rotations in the icy-plane. 

It is easy to see that one such block may be put on top of another in 

(26 - l)(2w - 1) + (b + w - l) 2 (1.1) 

different ways if b ^ w and in 

(26 -l) 2 (1.2) 

different ways if 6 = w. 

With w = 4 and 6 = 2 we get 46 possibilities, and note (as depicted in blue 
on Figure ^) that 2 of these are symmetric. Thus, letting H2xi(n,m) denote 
the number of ways to build a building of height m with n 2 x 4 LEGO blocks 
one then clearly has 

H 2x4 (n,n) = ±(46 n - 1 +2 n - 1 ) 

Note that i?2x4(6, 6) = 102981504, so that in fact LEGO's computation is off by 
four. 

By combining results of computer-aided enumerations with elementary com- 
binatorics one can further establish 

H 2 xi(n, n - 1) = 46 n ~ 4 (-89115 + 37065n) + 2™~ 4 (-8 + 16n) (1.3) 

for n > 3 and 

H 2X 4{n,n-2) = 2 n ~ 7 (1785 - 825??. + 256n) (1.4) 
+46™- 7 (-918674675- 5330182078n+ 1373814225n 2 ) 

for n > 5, but as the problem is rather hopelessly non-markovian there seems to 
be no way to give formulae for the number of buildings of relatively low height, 
or indeed for the total number T2 X i(n) of contiguous configurations, counted 
up to symmetry. Although symmetry arguments and other tricks can be used 
to prune the search trees somewhat, we are essentially left with the very time- 
consuming option of going through all possible configurations to determine these 
numbers, which even with efficient computers seems completely out of range for 
numbers such as T 2X 4(12). A sample of our results may be seen in FigureEl 

These results seem to indicate, as shown on Figure that T2y.i{n) grows 
exponentially in n. In this paper we will show that this is indeed the case, and 
give upper and lower bounds on the rate of growth - the entropy of the blocks. 
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Figure 2: i?2x4(^ : Tn) for m, n < 6 
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Figurc 3: Rate of growth of T 2 x4 with scmilogarithmic plot 

2 Entropy of b x w blocks 

It is the goal of the present section to prove that the following definition makes 
sense: 

Definition 2.1 The entropy of a b x w LEGO block is 

,. logCZbx^n)) 
Sbxw = hm (2.5) 

n — >oo Tl 

We let h bxw = exp(sb Xw ). 

That the limit exists is by no means obvious, except of course for s lx i = 0. 
We shall prove that this is the case in two steps, first establishing convergence 
in [0; oo] and then proving that the limit is finite. 

It is inconvenient and irrelevant for our theoretical considerations to identify 
buildings up to symmetry, so we establish definitcness in another fashion. Sup- 
pressing the block size from the notation, we will by A n denote all contiguous 
buildings containing [0;u>] x [0;6] x [0; 1] with the further property that there 
is no other block in M 2 x [0; 1] and no block at all in K 2 x [— 1;0]. Thus, the 
configuration can be thought of as sitting on a base block at a fixed position. 

We let a„ = #^4 n and note 
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Lemma 2.2 We have 

Um log(T bxw (n)) _ Um log(a„) 

n >oc fi n »oo fi 

in the sense that if one limit exists, so does the other. 
Proof: The claims follow immediately by the inequalities 

T b x w(n - 1) < a„ < 4T bxw (n) 

The leftmost inequality follows by mapping each equivalence class of configura- 
tions with n — 1 blocks to a representative placed on top of [0; w] x [0; b] x [0; 1] 
and noting that this map is injective. The rightmost follows by mapping each 
configuration to an equivalence class and noting that this map is at most 4 — 1 . 
□ 

We now get 

Proposition 2.3 log(a„)/n converges in [0; oo] as n — > oo. 

Proof: One sees that a n +m > a„a m by noting that an injective map from 
A n x A m to A m +n is defined by placing the base block of the element of A m 
somewhere on the top layer of the element of A n . 

Hence log(a„) is a supcradditive sequence, and log(a„)/n converges to 

sup log(a„)/rt. 

nGN 

□ 

To prove that this limit is finite, i.e. that a n grows no faster than exponen- 
tially, we describe a surjective map associating to each function 

S n : {1, . . . , 2bw(n - 2)} — ► {-few, -bw + 1, . . . , bw - 1, bw} (2.6) 

an element of 

An U {FAIL}. 

Clearly the number of such functions grows only exponentially in n. We 
shall subsequently look closer at which functions do indeed lead to buildings 
with n LEGO blocks, and give much better estimates for h^xw than the obvious 
(2bw + l) 2bw . 

With a fixed enumeration of the studs and holes of a ft x w LEGO block by 
the numbers 1, . . . , bw, a map of the form l|2.6|l gives rise to an element of A n , 
or the symbol FAIL, as follows. 

Take one LEGO block and call it block 1. Then read S n (l), . . . , S n (bw) from 
left to right to specify what to build on top of block 1 as follows. If > 0, 

take another LEGO block and place it parallely to block 1 with hole S n (l) on 
top of stud 1. If SVi(l) < 0, take a LEGO block and place it orthogonally, rotated 
+90°, to block 1 with hole —S n ({) on top of hole 1. In both cases, give the new 
block the number 2. If 5 n (l) = 0, do nothing. Then proceed to read 5„(2) to 
see what, if anything, to place on stud 2, and so on until S n (bw). Enumerate 
the blocks as they are introduced. 

Terminal state 2.4 If at any point a block collides with one which has already 
been placed, the procedure terminates with FAIL. 
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These steps will result in the placing of between and bw blocks on block 1. 



Terminal state 2.5 If at any point all n blocks have been placed, consider the 
unread values of S n . If they are all 0, the procedure terminates successfully 
with an element of A n . If not, the procedure terminates with FAIL. 

Terminal state 2.6 If, after reading the specihcations for the first m < n — 1 
blocks, no block m+1 has been introduced, the procedure terminates with FAIL. 

We may now assume that a block 2 has been introduced and look at S n (bw + 
1), . . . , S n (2bw) which will specify what to build on top of this block, if anything, 
in the same way that <SVi(l), . . . , S n (bw) specified what to build on block 1. A 
positive number at S n (bw + 1) will result in the placing of a block on stud 1 
of block 2 parallel to block 1, a negative number at S n (bw + 1) will result in 
the placing of a block on stud 1 of block 2 orthogonal to block 1, etc. We 
proceed in the same way for blocks 3, . . . , n — 1, but now read 2bw values where 
S n ((2m — 4)6u> + 1), . . . , S n ((2m — 3)bw) will specify what to put on top of block 
m, and S n ((2m — 3)bw + 1), . . . , S n ((2m — 2)bw) will specify what to put on 
underneath it in an analogous way. 

Terminal state 2.7 If at any point a second block is placed at the level M 2 x 
[0; 1], the procedure terminates with FAIL. 

Terminal state 2.8 If S n has been read to the end, consider the number of 
blocks placed. If it is less than n, the procedure terminates with FAIL. If not, 
it terminates successfully with an element of A n . 

We repeat this until one of the terminal states are reached. 

If the procedure does not fail, it will result in a building of n contiguous 
blocks, and clearly any such building may be constructed in this way. Thus, the 
number of possibilites for maps S n dominates a„, as desired, and we have: 



Theorem 2.9 The limit in lj2.5J) exists for any block dimension b 



x w. 



We can give general bounds of hbxw, but as wee shall see below in the case 
w = 4, b = 2, these bounds can in general be rather dramatically improved. 



Theorem 2.10 Ifb^w we have 

in - ~\\ 2 < h,.. < 

{2bw - 1)26x0-1 



(2b -l)(2w-l) + (b + w- If < h bxw < -A- >.-.. (2.7) 



Ifb = w we have 



( b 2\2b 2 + l 

(2b-lf<h bxw < {b ^ i)2b2 _ 1 (2.8) 



Proof: By (|1 . 1|> we clearly have 

a„ > ((26 - l)(2w - 1) + (b + w - l) 2 )"" 1 

when 6 ^ w, and this gives the lower bound in that case. For the upper bound, 
note that a function S n with nonzero entries at m locations will yield FAIL if 
m 7^ n — 1. Hence we get 
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Figure 4: Enumeration of studs and holes 



By Stirling's formula we then get 



< 



{2bw(n - 2))! 



(n - l)!((2btu - l)n - Abw + 1)! 
(2bwn)\(2bwT 
n\((2bw - l)n)\ [ ' 
{2bwn) 2bwn (2bw) r 



(2bwY 



v ' v ' 0(n- L/z ) 

n n ((2bw - l)n)( 2bw - 1 ) n 

{2bwf hwn (2bwY _ 1/2 

{2bw-l)( 2bw - 1 ) n [ ' 

(2bw) 2bw+1 



{2bw - 1 



|26tu-l 



0(n-^ 2 ) 



from which the claim follows directly. 

The square case follows similarly by l|1.2fl and by noting that functions S n 
may be chosen non-negative. □ 

Example 2.11 We enumerate the studs of a 2 x 4 LEGO block according to 
Figure^ Now consider functions {!,..., 32} i— » {—8, . . . , 8} given by 



(0, 5, 0, 0, -4, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, . . . , 0, 0, . . . , 0) 

12 3 



(0, -1, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, ... , 0, 0, 0, 0, -1, 0, 0, 0) 

12 3 



(0,. ..,0,0,. ..,0,0,... ,0,0,. ..,0) 

12 3 



(0, -1, 0, -1, 0, 0, 0, 0, 0, . . . , 0, 0, 0, 2, 0, 0, 0, 0, 0, ... , 0) 

1 2 3 

(0, -1, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, . . . , 0, 0, . . . , 0) 



(2.9) 
(2.10) 
(2.11) 



(1, 1, 0, 0, 0, 0, 0, 0, 0, . . . , 0, 0, . . . , 0, 0, . . . , 0) 

12 3 

(0, -1, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0,0, ...,0,0,0,0, -1,0,-1,0) (2.12) 

12 3 

(2.13) 
(2.14) 



(2.15) 



where all ellipses indicate six consecutive zeros. The functions H2.9(l and (|2.10(l 
give rise to the buildings depicted on Figure The remaining five functions 
give simple examples of functions resulting in the procedure failing at Terminal 
state rO \TR mi VTA and El respectively. 
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Figure 



5: Buildings associated to (|2.9[) and (|2.10|) 



3 Improved upper bounds 

In this section we shall describe methods to improve the upper bound on rifrxw 
given in Theorem l2.10l They apply to any dimension, but as they are somewhat 
ad hoc we shall concentrate on our favored dimension 2x4 and leave other cases 
to the reader. 

From Theorem l2~TUl wc know that fo x4 < 16 17 /15 15 < 647.02. We shall give 
a simple improved estimate leading to fo X 4 < 203.82 and a somewhat more 
complicated one leading to fo X 4 < 191.35. Besides being easier to state, the 
simpler estimate has applications in producing statistical estimates for a„ for 
relatively large n. 

Note that the surjective map associating buildings (or FAIL) to certain maps 



is very far from being injectivc. Wc have already employed the fact that unless 
the number of nonzero values is n — 1, the function is mapped to FAIL. But we 
may also use that the placement of a block onto another may be indicated in 
I different ways, where t is the number of studs of the lower block which are 
inserted into the upper block. Restricting attention to maps where placements 
are indicated in a fixed way will not affect surjectivity of the map. 

Any partition of the 46 positions in Figure into 8 sets Vi,...,Vs with 
the property that any position in Vi employs stud i of the lower block can be 
used to improve the upper bound. One uses the convention that a position in 
Vi is always indicated by a symbol at stud i, thus restricting the number of 
possibilities. 

Another restriction is available when specifying what to add to block m for 
to > 2. If we keep track of how block m was introduced, we know a priori 
that one hole or one stud of it has already been used, thus eliminating at least 
16 out of the 46 possibilities on the relevant side of the block. Dividing up the 
remaining 30 positions as above, we get 64 sets V\ with the property that V\ = 
and that V{ , . . . , V\ is a partition of the 30 positions which do not employ stud 



S n : {!,..., 16(n- 2)} 



{-8, -7,..., 7, 8} 



3- 



Theorem 3.1 We have 




(3.16) 



and, consequently, that fo x4 < 6 • 13 13 /12 12 < 203.82 



7 




Figure 6: Even distribution 

Proof: We partition the 46 positions into 8 sets, each consisting of 6 or 5 con- 
figurations, as indicated by the rows of FigurcEl Thus on indices related to one 
side of blocks 1, . . . ,n — 1, we need only allow for 6 different symbols. 

On the other side of blocks 3, . . . , n — 1, we can do even better, as outlined 
above. We leave to the reader to check that the 30 positions may be distributed 
evenly over 5 studs. Hence, H3.16J) is established, and the remaining claim follows 
by Stirling's formula as in Theorem 12. 101 □ 

It turns out - somewhat countcrintuitively? - that uneven distributions of 
the positions give slightly better estimates than what we obtained above. We 
have not carried out a systematic analysis and can not claim that the distribu- 
tion leading to Theorem 13.31 is optimal, but trial and error with the following 
proposition make us believe that there is only marginal room for improvement 
by this method. 

If (ai, . . . , etg) and (b\, . . . , bg,) are tuples of integers, we write (oi, . . . , as) < 
b%) if there is a permutation a of {1,...,8} with the property that 
Oo-(i) < bi for each i. 

The methods leading to the following result are surely known. 

Proposition 3.2 LetVi andV\ be partitions of the sets of positions as outlined 
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above, and assume that 

(#Vi,...,#V a ) < (oi,...,ob) 
(#^',...,#P|)<(6i,...,6 8 ), i€{l,...,8}. 

fo(y) = (y + 01) • ■ • {y + ob) P(y) = Po(y)(y + h)---(y + b s ) 

we have that a„ is dominated by the coefficient of y 15n ^ 31 in (Po(y)) 2 (P {y)) n ~ 3 
and that 

^2X4 < —35- 

where xq is the largest real root of 

Q{x) = l5P(x) - xP'(x) 

With the even distribution described above we get xq = 72, which, since 
P(72)/72 15 = 6 • 13 13 /12 12 is consistent with Theorem O Using that in fact 
. . . , #Vs) < (5,5,6,6,6,6,6,6) we may improve the estimate Oil fT>2x4: 
sligthly to 198.57. 

However, we can do even better with very uneven distributions: 

Theorem 3.3 h 2/A < 191.35 

Proof: There exist partitions with 

(#7>i,...,#P 8 ) < (16,15,7,5,2,1,0,0) 

< (15,7,4,3,1,0,0,0), je{l,...,8}. 

so by Proposition 13. 21 we are lead to consider 

P(x) = -x 5 (x + I5)(x + 7){x + l)R(x) 

where R(x) is the polynomial 

.T 8 -23x 7 -2056x 6 -38700x 5 -332657a; 4 -1504645a; 3 -3645736a; 2 -4392600a;-2016000 

which has a largest real root which is approximately 65.05. The estimate follows 
by Proposition ^. 21 □ 

4 Improved lower bounds 

It follows from Theorem 12 . 1 01 that /i2x4 > 46. We shall in this section improve 
this estimate to /i2 X 4 > 78.32. 

We let B n C A n +i denote the set of LEGO configurations as above consisting 
of n + 1 blocks and such that both the top and the bottom layer consists of a 
single block. Setting b„ = #£>„ we then clearly have 

a n < b„ < a„ + i 



9 



Figure 7: Two bottlenecks 



and hence, as in Lemma 12 . 21 



^2x4 = ex P li m 




log b 



n 



n 



) 



(4.17) 



We say that a configuration c in B n has a bottleneck at height z G N if c 
has exactly one block in the layer K 2 x [z; z + 1]. By convention the top and 
bottom blocks are not bottlenecks. This ensures that removal of a bottleneck 
decomposes c into two configurations c' and Cq one of which, say c , contains the 
bottom block of c. Re-inserting the removed block in c' yields a configuration d 
in some B n with the inserted block as the top block. Re-inserting the removed 
block into c ' yields, after a translation, a configuration c" in B n - m with the 
inserted block (translated) as the bottom block. Evidently, we can reconstruct 
c in a unique fashion from (c', c"). Repeating this decomposition procedure we 
conclude that any configuration in B n with exactly k > bottlenecks can in a 
unique way be decomposed into a sequence (c^\ . . . ,d k+1 ^) of configurations 
such that c'*) € S mi has no bottlenecks and mi + • • • + nik+i = n. 

Letting C n denote the subset of B n consisting of configurations without bot- 
tlenecks we obtain in this way a one-to-one correspondence between elements of 
B n and those of 



oo 




miH \-rn k+1 =n 




c. 



x • • • x C, 



(4.18) 



Let now 



and let ^ and V'o denote the generating functions 



oo 



oo 





(4.19) 
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It follows from l|4.18[l that 



^)=f> (z))< = -^°^. (4.20) 

From the definition of ip(z) and h 2 xi it follows that ip is analytic in the disc 

D = {z\\z\<{h 2 ^)- 1 } 

and, since b n > 0, that tp is non-analytic at z = {h 2 y.i)~ 1 - From H4.20JI we hence 
conclude that 

|VoO)| < 1 for |z| < (fox*) -1 (4.21) 

In particular, we get 

Cl(/l 2 x4) _1 + • • • + C n (fo x4 )~" < 1 

which gives our claimed lower bounds on fo X 4, depending on the number of 
terms n on the lefthand side. 

We shall describe in detail how to get the first order of improvement of the 
estimate from Theorem 12 .101 As evidently C2 = we turn to C3 for this. 

The configurations contributing to C3 have one bottom block, one top block, 
and two blocks in between. The number of ways of placing two blocks on top of 
the bottom block is rather easily seen to be 480, so the number of configurations 
where the two middle blocks are both attached to the bottom block is 2-46-480 — 
4730, where 4730 is the number of configurations where the middle blocks are 
both attached to the top block as well as to the bottom block. The remaining 
configurations are those where the middle blocks are both attached to the top 
block but only one of them to the bottom. This number is seen to be 2 • 46 • 
480 - 2 ■ 4730. Thus we have 

c 3 = 4 • 46 • 480 - 3 • 4730 = 74130 

and P3{h 2xi ) — 1 where 

P 3 (x) = 46a; + 74130a; 3 

This gives 

fo x4 > 64.06 
which can be improved as follows. 

Theorem 4.1 fo x4 > 78.32 

Proof: Computer-aided computations give 

c 4 = 867346 
c 5 = 318434429 
c 6 = 18335373238 

so we have that -P6(^2x4) — 1 where 

P 6 {x) = 46x + 74130a: 3 + 867346a; 4 + 31843442&E 5 + 18335373238a; 6 
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which gives Zi2x4 > 76.67. 

To improve the estimate we prove that 

c„ +2 > 1248c„ (4.22) 

for n > 6. To see this, we devise 1248 different ways to construct an eiement c' 
of C„+2 from an element c of C n , in such a way that the original configuration 
can be recovered from the resulting one. 

Of these 1248 configurations, 480 are gotten from a fixed configuration d of 
two blocks sitting on one base block b by identifying the base block of c with the 
block at the second level of d which meets the stud of lowest index on b according 
to the enumeration of Figure 0] We get the configuration c' by rotating 90°, 
if necessary. The remaining 768 configurations are gotten by placing one block 
underneath the base block of c, and placing one more block at the level of this 
original base block, such that these two added blocks do not meet. A computer 
search shows that there is always at least this number of ways to do so since 
there are at least two blocks at level 1 of c but only one at level 0. We get the 
configuration c! by translating the configuration upwards and rotating 90°, if 
necessary. 

To reconstruct c from c', one first sees how many blocks are attached to the 
base block of c'. If there are two, c is gotten by discarding the base block of c' 
and the block at the next level sitting at the highest index of it, translating down 
and rotating 270°, if necessary. If there is only one, c is gotten by discarding 
the base block of c' and the block at the next level which does not meet that 
block, translating down and rotating 270°, if necessary. 

By repeated application of (|4.22|) we get c 6+2 fc > 1248 fc C6 and 05+2/; > 
1248 fc c 5 , so that r(h 2 xi) < 1 with 



r(x) = P 6 (x) 



C 5 X 7 + CqX 8 



1/1248 -a; 2 ' 

leading to the stated lower bound. □ 

5 Concluding remarks 

We do not at present have the software nor the computer power to perform 
numerical experiments to get a good idea of the true value of /i2x4- Our best 
guess, based mainly on data achieved by Abrahamsen on the presumably closely 
related case of 1 x 2-blocks would be that the number is rather close to 100. 
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